Discrete solitons in coupled active lasing cavities 
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We examine the existence and stability of discrete spatial solitons in coupled nonlinear lasing cavities 
(waveguide resonators), addressing the case of active media, where the gain exceeds damping in the linear 
limit. A zoo of stable localized structures is found and classified: these are bright and grey cavity solitons 
with different symmetry. It is shown that several new types of solitons with a nontrivial intensity dis- 
tribution pattern can emerge in the coupled cavities due to the stability of a periodic extended state. The 
latter can be stable even when a bist ability of homogenous states is absent. © 2012 Optical Society of America 
OCIS codes: 190.1450, 190.4390, 190.4420 



Cavity solitons have attracted appreciable attention as 
their stability and high robustness make them excellent 
candidates for multipass optical delay lines, spatial in- 
formation encoding, and manipulation/reconfiguration 
based on the transverse distribution of light [1,2]. These 
objects represent a particular realization of the dissipa- 
tive soliton notion - localized solutions of nonlinear equa- 
tions existing due to the balance between attenuation, 
gain, nonlinearity and diffraction [2] . Recently there has 
arisen a significant interest in studying discrete cavity 
solitons (DCS) [3-8], which are localized excitations in 
coupled arrays of nonhnear cavities. Notably, the exis- 
tence and stability properties of DCS can be essentially 
different from those of the continuous counterparts, and 
moreover, the continuum limit of such solutions can be 
absent at all [3-5] . At the same time the self-organization 
and emergence of localized structures in lasing cavities 
are also the subjects of intensive theoretical [2,8] and ex- 
perimental [9,10] study. The important difference of las- 
ing systems is that there one can have solitary solutions 
maintained by a simultaneous balance between both lin- 
ear and saturable gain/loss impact [2,6,8]. In this Letter 
we consider DCS in the system of coupled lasing cav- 
ities, taking into account linear loss, Kerr nonlinearity 
and gain saturation. For the first time to our knowledge, 
we address the localized excitations in discrete active 
systems, where in a low amplitude regime the gain term 
dominates over the attenuation. In this Letter we show 
that, due to the presence of external driving, this system 
can support versatile families of stable localized solutions 
unless the gain is bigger than a particular value. We also 
present a new type of DCS^ which has not been studied 
so far in other discrete coupled-cavity systems. 

We consider the array of waveguides and will use the 
model developed in [4] taking into account the interac- 
tion between the neighboring cavities only. Characteriz- 
ing the field in each cavity by slow varying amplitude An 
we write down the resulting normalized equation for the 



field amplitude inside the n-th optical resonator as: 

+ C(^„+i+A„_i-2A„) =P-e^'=". (1) 

In Eq.(l) C characterizes the strength of coupling, A = 
A — 2C is the detuning of the pump frequency from the 
resonant one (we explicitly extracted the term 2(7 to 
have a second difference); parameter ^ > describes the 
linear losses, P > is the amplitude of the driving field 
(normal incidence is considered, = 0), a characterizes 
the strength of the nonlinear (Kerr) polarization, and 
7 > is the lasing gain. We are engaged with the case of 
active lasing media^ and thus inequality j > S holds; this 
implies that in the absence of driving field (P = 0) the 
trivial solution = is modulationally unstable [8]. 
The field An is normalized to the lasing saturation inten- 
sity Iq: Anl^fh) the time t is also normalized by 
the cavity relaxation time (for the exact normalization 
in Eq.(l) see [4,7]). For presenting the results we opt for 
the following typical values of normalized input param- 
eters: A = 3.5, (5 = 0.4, 7 = 1.7. Further we focus on the 
discrete solutions that do not have continuum counter- 
parts [4,7], and deal with a weak defocusing nonlinear- 
ity: a = —0.2. So, the solutions can exhibit sharp spatial 
changes with scales of the order of the inter-cavity sep- 
aration. We use (unless otherwise stated) the relatively 
weak coupling strength C = 0.15. We notice that by ap- 
plying the symmetry transformations [7] one can readily 
compose the counterparting solutions for other values of 
parameters starting from the families given below. 

First we address the extended static states of Eq.(l). 
For the homogeneous (H) pattern, setting An = A, one 
finds the so-called response curve P(|A|): it can have ei- 
ther one or three solutions depending on the system pa- 
rameter values, and in the latter case one can have bista- 
bility between the low-intensity (HI) and high-intensity 
(112) homogeneous states. The curve P(|A|) for the set 
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Fig. 1. (Color online.) Response curves for homogeneous 
and periodic solutions with stable regions marked by 
thick lines. Pane (a): The curve P(|A|) for the homo- 
geneous H-solutions, with stable lower (HI) and higher 
(H2) states (the dash-marked top of HI is also unsta- 
ble for C > 0); below this curve we show the spectrum 
of Re [A] (shadowed region), the boundaries (7 = and 
q = 71 are drawn by sohd lines. Pane (b): The curves 
P{\^max\) and P{\Amin\) for the periodic inhomoge- 
neous I-solution, the direction of 'evolution' from the 
curve coalescence points (marked by ||) is identified with 
arrows. The inset shows the I-solution profile correspond- 
ing to the point on the snakes. The horizontal dash- 
dot lines and the highlighted parts on the P-axis indicate 
the bist ability region for homogeneous solutions. 



of parameters chosen is given in Fig. 1(a) and it is mul- 
tivalued for some diapason of P. 

One can also find the inhomogeneous periodic static 
solutions of Eq.(l) (we label these as I-solutions), hav- 
ing a 2-site period: we set An, etc. = Amax^ and 

An^l, 7 etc. = Amin-) Amax 7^ Amiri') and then SOlvC 

the system for these two complex fields. The results 
are given in Fig. 1(b) in the form of two (equivalent) 
snaking response curves: P{\Amax\) and P{\Amin\)' The 
I-solutions are discrete and do not exist in the continuum 
limit. In general, infinitely many extended solutions with 
more complex spatial patterning can be found but here 
we restrict ourselves to these three simplest families. 

We are aiming to study solitons forming on H or I 
backgrounds and so the stability of these states is of 
importance. To analyze the stability we linearize Eq.(l) 
and solve the eigenvalue problem for small excitation 
an{t) = a(i)e'^^+*^^ + a(^2)^^*^~^^^ on the background A. 
The background is stable if Re [A] < for all A. 

We see. Fig. 1(a), that for the parameters chosen we 



have a region of bistability (highlighted). Note, that the 
modulational instability of the part of Hl-branch for low 
P and \A\ values is a characteristic feature of active cav- 
ities, and is absent for passive case [4,6,7]: the stability 
region grows towards lower values of P as the 'activeness' 
of the medium diminishes, i.e. when 7 tends to 5. For any 
5 there is a threshold value of gain, 7max(^^, ^1(^)1 where 
the the only stable H-state is the high-intensity one, H2. 

The stable region of I-state corresponds to the high- 
lighted snake parts in Fig. 2(b). The stability regions of 
HI and H2 partially overlap with that of the I-solution 
but the I-state can he stable when HI is unstable, i.e. 
in the region where one does not have the bistability 
of H-states. We emphasize that this nonlinear extended 
I-state has not been studied before for the systems of 
coupled optical cavities [3,4,7], and its existence and 
stability allow composing new stable DCS families. 

Static DCS can now be found in the form of (ho- 
moclinic) connections between the aforementioned three 
stable states: homogeneous ones HI and H2, Fig. 1(a), 
and a periodic inhomogeneous I-state, Fig. 1(b). In Figs. 2 
and 3 we protocol our DCS found in the model Eq.(l), 
presenting these on the snaking diagrams: each point on 
a snake corresponds to a particular DCS. The DCS can 
be divided into two subclasses according to whether the 
intensity in the soliton center is higher (bright solitons. 
Fig. 2) or lower (grey solitons. Fig. 3) compared to the 
background [4, 7]. It is convenient to characterize the 
bright/grey DCS by the maximal/minimal amplitude of 
the soliton. The simplest DCS correspond to the con- 
nection between two homogeneous states: H1-H2-H1 for 
bright DCS, i.e. 'humps' of intensity residing on the HI 
background (see Figs. 2 (a)- (c) and the profiles of stable 
solution at panes corresponding to the points • '1', '2' 
and '3'), and H2-H1-H2 for grey DCS, i.e. an intensity 
'dent' in the H2 background (not shown). Both these 
types of DCS can be stable in the active media described 
by Eq.(l), see the highlighted regions on the snaking 
curves in Figs.2(a)-(c) for the stable bright DCS. The 
solitons can exist only up to a finite value of coupling: 
in Fig. 2 (a) the C(\Amax\) dependencies for bright DCS 
turn down after passing a maximum point, and the same 
situation holds for the grey solitons, i.e. these solutions 
do not exist in the continuum limit C ^ 00. The snaking 
diagram in terms of P{\Ajnax\) for the H1-H2-H1 bright 
DCS is given in Fig. 2(b), but since the stable winds of the 
PdA^aa^D-snake are hardly distinguishable for this soli- 
tonic family, a more representative diagram can be plot- 
ted in terms of the solution mass, M = \An — ^bP, 
with Ab being the background field value. On this plot 
[Fig. 2(c)], the hierarchy with respect to the growth of 
the number of excited DCS sites is clearly seen. 

The new families of DCS correspond to connections 
of the homogeneous states with the I-state: bright Hl- 
I-Hl, and grey H2-I-H2; see Figs.2(d), (e) and Fig.3, 
correspondingly, for the snaking diagrams in terms of 
P{\A 

max I), M{\A 

max 

I) and P{\Amin\)' As one can see, 
these DCS have a nonuniform central excited part cor- 
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Fig. 2. (Color online.) Snaking diagrams (stable regions 
are highlighted) and solution profiles for bright DCS. (a) 
The curve Cd^maxl) for fixed P = 4.47 corresponding 
to the 1- and 2-site H1-H2-H1 solutions. Panes (b) and 
(c) are the diagrams for the site-centered H1-H2-H1 so- 
lutions in terms of P{\Amax\) and M(\Amax\)-, panes (d) 
and (e) - the same for the Hl-I-Hl solutions. The small 
panes right to panel (a) and below the figure show the 
particular stable distributions referring to the points on 
the snakes marked as • 'l'-'5'. 



responding to a fragment of the I-state embedded into 
either HI or H2 background (see the profiles for points 
• '4', '5', and ■ '6'). These 'nontrivial' solitons are sta- 
ble for some range of parameters, see the highlighted 
parts of the snakes. Note, that grey DCS of this type 
can be stable in the region, where solitons referring to 
the connections of the H-state are unstable. Fig. 3. It is 
also possible to compose yet another new type of stable 
DCS in the form of a uniform HI or H2 segment embed- 
ded into a periodic inhomogeneous I-background, due to 
the intersection of the stability regions for the H- and 
I-states. 

To conclude, we have found a zoo of localized static 
solutions corresponding to stable DCS in coupled las- 
ing cavities; bright DCS are protocolled in Fig. 2 and 




Fig. 3. (Color online.) Snaking diagram P{\Amin\) for 
the grey DCS corresponding to the H2-I-H2 connection. 
The inset shows the particular distribution referring to 
the point marked as ■ '6'. 



grey ones in Fig. 3. Aside from 'conventional' DCS, cor- 
responding to the connections between homogeneous 
states HI and H2, we have found a new type of DCS 
involving a periodic inhomogeneous I-state, which has 
not been found before in optical cavities. The existence 
of the great variety of stable DCS paves the way to the 
more versatile and sophisticated patterning of transverse 
light distribution. At the same time our study poses an 
interesting question of whether the 'nontrivial' DCS, like 
those exemplified by solutions • '4', '5' and ■ '6', can be 
stable in other systems of coupled optical cavities , like 
a discrete Lugiato-Lefever model [4] . Notably, these soli- 
tons can be stable when the bistability of H-states is 
absent. 
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